wage_and_price_file = '../../../data/generated/matlab/OptimTax/gpops/optimum_vars.csv';
knitro_path = 'PATH_TO_KNITRO';

%% compute optimal robot tax, keeping tax progressivity fixed

% get the initial tax progressivity from here
init_file = '../../../data/generated/matlab/OptimTax/hsv/20191107_nochoice_adjust_tau_hsv.csv';

% wage_and_price_file is needed to get the price levels?

results = optimtax_hsv('adjust_robot_tax',knitro_path, 'wage_and_price_file', ...
             wage_and_price_file,'compensation_init_file',init_file);
         
         
%% use results from step above and compute compensation
init_file =   ['../../../data/generated/matlab/OptimTax/hsv/' ...
               '20191107_nochoice_adjust_robot_tax.csv'];
welfare_targets_file = init_file;

results_comp = optimtax_hsv('baseline',knitro_path,'compute_compensation',true,'welfare_targets_file',welfare_targets_file,'suffix','optimal_robot_tax_hsv','compensation_init_file',init_file);